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Abstract. We propose a mathematical theory for the refocusing properties observed in 
time-reversal experiraents, where classical waves propagate through a medium, 
are recorded in time, then time-reversed and sent back into the medium. The 
salient feature of such experiments is that the refocusing quality of the time- 
reversed reemitted signals is greatly enhanced when the underlying medium 
is heterogeneous. Based on the Wigner transform formalism, we show that 
random media indeed greatly improve refocusing. We analyze two different 
types of random media, where in the limit of high frequencies, the Wigner 
transform satisfies a random Liouville equation or a linear transport equation. 



Renversement du temps pour la propagation d'ondes clas- 
siques en milieu aleatoire 

Resume. Dans des experiences de renversement temporel, oii un signal acoustique 
qui se propage dans un milieu sous-jacent, est enregistre en temps, puis 
subit un renversement temporel et est reemis dans le milieu, il a ete 
observe que le signal refocalise d'autant mieux a 1' emplacement de la 
source initiale que le milieu est heterogene. Nous donnons une expli- 
cation mathematique de ce phenomene fondee sur le formalisme de la 
transformee de Wigner. Nous obtenons deux exemples de refocalisation, 
en fonction du type de milieu aleatoire considere, selon que la trans- 
formee de Wigner, a la limite des hautes frequences, satisfait soit une 
equation de Liouville avec coefficients aleatoires, on une equation de 
transport. 



Version frangaise abregee 

Dans cette note, nous nous interessons a la propagation d'ondes classiques en milieu aleatoire et a 
leurs proprietes de refocalisation apres renversement temporel. Dans les experiences conduites par 
M. Fink Q et ses collaborateurs, une source localisee emet des ondes acoustiques a tr avers un milieu 
heterogene. Le signal est enregistre en temps par des capteurs localises en espace et positionnes loin 



de la source, puis renverse en temps et enfin retransmis dans le milieu (ce qui est arrive en dernier 
repart en premier). 

Dans ces experiences, il a ete observe le phenomene suivant, qui a premiere vue pourrait paraitre 
contre-intuitif : le signal reemis refocalise a I'emplacement de la source initiale d'autant mieux que le 
milieu est heterogene, alors que la qualite de la refocalisation est tres mauvaise en milieu homogene. 
L'explication reside dans la multiplicite des "chemins" menant de la source aux capteurs. Plus 
le milieu est aleatoire, plus il donne d'informations sur la source aux capteurs, pourvu que les 
enregistrements soient sufRsamment longs. 

Plusieurs travaux sur ce phenomene existent dans la litterature physique |y, y, |7[ • Les resultats 
mathematiques existants portent sur le cas mono-dimensionnel [^ et le regime de faisceaux etroits 
Ig]. Dans cette note, nous proposons une theorie qui explique le phenomene de refocalisation pour 
des ondes classiques a hautes frequences en milieu multi-dimensionnel. Notre approche, s'inspirant 
de ||2|, se fonde sur la transformee de Wigner, definie en @, qui permet I'analyse des correlations de 
noyaux de Green en des points voisins. Nous pouvons ecrire le signal retransmis apres renversement 
du temps en fonction de cette transformee de Wigner @. 

Selon les caracteristiques du milieu aleatoire, la transformee de Wigner satisfait differentes 
equations limites quand la frequence caracteristique des ondes tend vers I'infini. Nous analysons 
les deux cas oil nous obtenons a la limite soit une equation de Liouville (limite semiclassique) , soit 
une equation de transport de type Boltzmann (limite de couplage faible). Dans le premier cas, nous 
sommes capables d'obtenir une derivation rigoureuse mathematiquement. La formule de refocali- 
sation est donnee en (p|)-(p!2[). Dans le second cas, la derivation de la formule de reconstruction 
(H)"(0) 6st formelle en certains endroits. La theorie permettant d'obtenir I'equation de transport 
limite n'est complete que pour I'equation (scalaire) de Schrodinger |5|, Q. De plus elle necessite de 
moyenner sur les realisations du milieu aleatoire, bien que les experiences physiques montrent que 
ce soit super flu. 



1. Refocusing in time-reversal 

Time-reversal experiments g] can be briefly described as follows. A point source signal is sent 
through an inhomogeneous medium and is recorded by a spatially localized array of receivers- 
transducers. The signal is subsequently reversed in time and reemitted back into the medium; that 
is, the part of the signal that was recorded last is sent back first. The refocusing of the reemitted 
signal is then observed on the spot of the original source. 

Refocusing in a homogeneous medium is poor when only a few receivers are used to record the 
signal. The most striking and somewhat counterintuitive observation is that inhomogeneities in the 
medium enhance the refocusing effects. This is because waves recorded at the time-reversal array 
have sampled a larger part of the medium than in the homogeneous case and carry more information 
about the source location. Furthermore, the refocused signal in a random medium is self-averaging, 
that is, independent of the realization of the random medium. 

Many time-reversal experiments have been conducted in the recent past and studied theoretically 
], 0]. Mathematical studies have concentrated so far on the one space dimension case [^ and on 



the narrow beam approximation |g, 12|. In this note we present an analysis of the enhanced refocusing 
property in genuinely three dimensional heterogeneous media in the regimes of random geometrical 
acoustics and radiative transport. Our main tool is the asymptotic analysis of the equation for the 



Wigner transform |Q, |rT|. It turns out that in the high frequency hmit the refocused signal may be 
written as a convolution with a kernel F expressible in terms of the limit Wigner distribution. More 
precisely, this kernel is expressed @ in terms of the Fourier transform in k of the energy densities 
a(i,x, k) in the phase space @. The latter are smooth in k under evolution in random media but 
not in a homogeneous domain. Therefore the spatial localization of the kernel F is improved in a 
random medium compared to the homogeneous case. In the geometric acoustics limit, the Wigner 
transform asymptotically satisfies a Liouville or a Fokker-Plank equation similar to the one used to 
describe the refocusing of narrow beam signals |2|, ^]. In the radiative transfer regime, the Wigner 
distribution satisfies a linear Boltzmann equation. The rigorous passage to this limit is not complete 
for classical waves P, 0, 13, 14] and our results in this part remain formal. 



2. Reemitted signal and the Wigner transform 

The propagation of many classical waves is described by systems of m first-order differential equations 
of the form (repeated indices are summed over throughout the paper) 

^(x)^^+D^-^^ = 0, xGM^<i = 2or3 (1) 

at oxj 

with initial conditions u(0,x) = uo(x). In this paper, we address time-reversal for acoustic waves 
with d = 3 to fix notation. The corresponding theory for other types of waves will be presented 
elsewhere [j|]. In this context, u is the 4- vector composed of the three components of the velocity 
field V and the scalar pressure field p. The matrix A{x) = Diag(/9(x),p(x),p(x), k(x)), where 
p{x.) is the density of the underlying medium and k(x) its compressibility. The 4x4 matrices 
D^ have entries D^, = 6mASij + dudmj- The solution of (|) is u(t,x) = La G(t,x;y)uo(y)(iy. Here 
G{t, x; y) is the matrix- valued Green's propagator from y to x. It solves (||) with the initial condition 
G(0, x; y ) = 16{x — y ) where / is the 4x4 identity matrix. Our time reversal setting is as follows. At 
time t the wave field is truncated in a region fi, where the transducers-receivers are located. Then 
the direction of the acoustic velocity field is reversed and the signal is re-emitted into the medium. 
The back-propagated signal emanating from the source term uq at time can be written using the 
Green's propagator as 

u«(x)= / TG{t,x;y)rG{ty-z)xniy)xn{y')f{y-y')no{z)dydy'dz. (2) 

Here T is a matrix that models the time reversal process. It is given by F = Diag(— 1, —1, —1, 1), so 
that the velocity field v is replaced by — v and the pressure field p remains unchanged. The function 
Xn{y) niay equal 1 on fi and outside, or could be a more general function modeling amplification of 
the time-reversed signal by various receivers. We also allow for some blurring of the recorded signal 
before it is transmitted back. The convolution with a filter / describes this blurring. Introducing 
the adjoint Green's matrix G*, solution of 

^^i|p^ + ^{G.{t,x;y))D^A-\x) = 0, (3) 

with initial condition G*(0,x;y) = 6{x — y)A~^{x), we observe that FG^=(t,x;y)^(x)F = G(t,y;x). 

We now rescale our problem. Refocusing is expected to be important in the close vicinity of 

the source location. We take the initial source as a localized function uo(x) = S ( ^"^^° ) of finite 



amplitude and rescale the filter accordingly: — /( ^~^ )■ Here e ^ 1 is a small parameter that 
measures the ratio of the width of the initial pulse to the propagation distance L between the source 
location and the array of receivers-transducers. An observation point x is close to xq and we write 
it as X = xq + e$,. After a change of variables, equation (|2|) becomes 

uf (^;xo) = y" rG(t,xo + e^;y)G.(t,xo + ez;y')A(xo +ez)rS(z)xn(y)Xf7(y')/(^^^)dydy'ciz. 

IR3d 

The Wigner transform is a natural tool in the analysis of the correlation function of wave fields at 



Xn{y)xniy')fi^^^)dydy'. (4) 



neighboring points ||8|, [11|, |13[. In our context, we define it as 

We{t,^,k) = [ [/■e^k.z^(t,x-^;y)G.(t,x+^;y')-^ 
This allows us to recast the expression for the retransmitted signal as follows 

uf(^;xo)=/ TWs{t,xo + e^,k)e-'''<'~^^A{xo + ez)TS{z)dzdk. (5) 

yR2d z 

Let us define the space A{M.'^'^) as the subset of 5'(M^'*) of matrix-valued distributions ?7(x, k) such 
that /jgd supx ||?7(x, y)||(iy is bounded. Here fj{x, y) is the inverse Fourier transform of ?7(x, k) in the 
second variable only. We denote by A' the dual space to A. Then we have 

Lemma 1. The Wigner transform VFe(t,x, k) is bounded in C^{0,T;A'{M.'^'^)) independent of e 
provided that /(k) E L^{]R. ). As a consequence, it converges weakly along a subsequence e^ ^ to 
a distribution Ty(t,x,k) G C^{0,T;A'{M.'^'^)). 

The proof of this lemma is obtained by rewriting W^ as 

Weit,^,k)= f /(q)/ e*k.z^,(t,x-^;q)a.(t,x+^;q)-^dq, 

where G'£(t,x;q) = fj^aG{t,x;y)xn{y)e^^'^'^^^dy solves (||) with initial data xn(x)e^*'i'^/^I, and 
G'e*(t,x;q) = f^aG,{t,x;y)xn{y)e'''"'^'dy solves (|) with initial data xn(x)e*'i''/"^-i(x). The 
functions Gg and G^* are uniformly (in e and q) bounded in L^, and hence P, |l|] f^d e*'''^G(t,x — 
f ;q)G'*(t,x-h f ;q)(^ is uniformly bounded in C°{0,T;A'{R'^'^)). This implies the result of the 
Lemma. 

The Wigner distribution at time t = is given by 1^(0, x,k) = \xn{'^)\'^f{k.)A^^{-x.). The dispersion 
matrix L(x, k) = A~^{x)kjD^ has a double eigenvalue coq = that corresponds to vortical modes, 
and simple eigenvalues ti;i^2 = ic(x)|k|, c(x) = 1/\/p(x)k(x). As was shown in ^, |^] the limit 
Wigner distribution may be decomposed as 

2 

W{t, X, k) = ^ a%{t, X, k)b°bf + ai(t, x, k)b^b^* + a2(t, x, k)h^h'^* . (6) 

Here b are eigenvectors of the matrix L. If the source S(z) has no vortical waves, that is, if S(k) is 
parallel to k, then the limit of the back-propagated signal is 

u^(^;xo) = V / a™(t,xo,k)rb"(xo,k)b"^(xo,k)e-^''(^-«)^(xo)rS(z)dzdk (7) 



m=l 



[ F{t, i-z; xo)S(z)dz = (F(t, •; xq) * S)(0, 

J'Rd 



(8) 



with ai(0,x, k) = 02(0, x,k) = |xf7(x)P/(k). The quahty of the refocusing of the back-propagated 
signal is determined by the decay properties in ^ of the kernel 

F(t,^;xo)=V/ a„(t,xo,k)rb™(xo,k)b'"(xo,k)e^l'-^A(xo)r(ik. (9) 

In the limit where f = 6 and fi = M'^, we verify that u^(^;xo) = S(^). When all the information is 
propagated back, the refocusing is perfect, as is physically expected. However, when the medium is 
homogeneous and c(x) = cq, the amplitudes 01^2(^1 xo,k) = |xa(xo T cokt)p/(k) become more and 
more singular in k as time grows (their gradient in k grows linearly in time). The corresponding 
kernel F = Ffj given by (0) being a Fourier transform in k of a function with growing gradients 
decays more and more slowly (in ^) as time grows. This leads to poorer and poorer refocusing. 

3. Refocusing via the Liouville and Fokker-Plank equations 

In this section, we analyze the random geometric acoustics limit, where the correlation length 6 
of the random medium is small compared to the propagation distance L but large relative to the 
wavelength e. The mode a = 0} satisfies the Liouville equation 

dta^ + Vkc^5 • V^a^ - V^ug ■ V^a^ = 0, (10) 

where a;(x, k) = (cq + Voci(x/(5))|k|. The Liouville equation is obtained thanks to the a priori 
estimate of Lemma 1. We assume that ci(x) is a mean-zero stationary random process with a 
smooth rapidly decaying correlation function i?(x). In the limit (5 ^ we have Eja | -^ a(t,x, k), 
where a solves the Fokker-Plank equation nC 



a=i«-(*°-W«-j-«.K-v«<. (11) 

with the initial data a(0,x,k) = |xf7(x)p/(k). Here, k = |k|, and the diffusion matrix is 

a2i 



_-| />oo ff'^ f} 



-00 ^■^j^^m 

with k = k/fc. The operator on the right side of ( [Tl| ) is degenerate because Dij(\i.)kj = 0. However, 
it is hypoelliptic on X = R^ ^ ^k ^^^ ^^^^ possesses a smooth Green's function GL(t, x, k; x', k') on 
X so that 

a(t,x,k)= / GLit,^,k;y,k')\xniy)\^f{^',k')dydk'. 

Hence we can introduce the new filter Fl = F^ + Ff^ with 

Fl(t,^;xo)= / e^''Mi,xo,k)rbi(k)bi*(k)^ordk, (12) 

and Ff^ defined similarly. In the limit e — > and (5 — > we obtain the reconstruction (g[) for E {u } 
with F = Fl- Now the function a(t,x, k) is smoother in k than in the homogeneous case and thus 
the filter Fl decays faster than Fh- This produces sharper refocusing and eliminates spurious Fresnel 
zones. The process a {t, x, k) does not converge to the deterministic function a{t, x, k) in probability 
pointwise in x and k. However, for any test function (/)(x, k) G 5(]R^ ) the process </> (t) = {a, ,(/)) 
converges to (a, cj)) in probability. That means that the refocused signal u^ is deterministic. This 
self-averaging of the refocused pulse follows from slight modifications of the analysis in |1C, 12 1. 



4. Refocusing via the transport and diffusion equations 

The radiative transport regime arises when the random medium has weak fluctuations at the wave- 
length scale. This case is more complicated than the one treated in Section 3 because the waves 
now fully interact with the medium. Propagation of wave correlations is described in terms of a 
Boltzmann-type radiative transport equation ||5|, |l^, |14|. The random fluctuations of the density 
and compressibility are assumed to be of the form /Je(x) = po + V^Piij)i ^£(x) = kq + \/eKi(j), 
where po, kq are constants, and pi(x), ki(x) are mean-zero stationary random processes. A rigorous 
derivation of the transport equation for the Wigner distribution is only available for Schrodinger 
equations |^, |l^, Based on the formal analysis in |l^ for hyperbolic systems, we can still conjecture 
that ai^2(i)X, k) solve a linear transport equation, at least after ensemble averaging over realizations 
of the random fluctuations. It is moreover known that for long distances of propagation and large 
times, the transport solution becomes independent of the direction k and is approximated by that 
of a diffusion equation: 

dta{t, X, k) - D{k)A^a{t, x, k) = 0, (13) 

where D{k) is a diffusion coefficient that depends on the power spectrum of the random fluctuations 
(see [^] for instance). When the filter / is radially symmetrical, the solution of ( p!3| ) may be 
written as 01^2(^1 x. A;) = ijj{t,'x.,k)f{k). Here 'ilj{t,:x.,k) is the solution of ( [l3|) with initial data 
■0(0, X, A;) = |xr2(x)p. We then obtain that the new filter Fb is scalar and is given by 

FB(t,C;xo)=/ e^^-^V'(t,xo,A:)/>)dk. (14) 

The retransmitted field (in the weak sense and after ensemble averaging) is given now by (|8|) with 
F = FbF Note that similarly to the random geometrical optics case treated in Section 3 the 
amplitude a(i,x, k) is smoother in k than in the homogeneous case and the filter Fb decays faster 
than Fh, thus produces a sharper refocusing than in a homogeneous medium. 
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